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ABSTRACT 

We have begun a program to discover high-redshift supernovae (z ~ 
0.25-0.5), and study them with follow-up photometry and spectroscopy. We 
report here our first discovery, a supernova at z = 0.458. The photometry 
for this supernova closely matches the lightcurve calculated for this redshift 
from the template of well-observed nearby Type la supernovae. We discuss 
the measurement of the deceleration parameter q using such high-redshift 
supernovae, and give the best fit value assuming this one supernova is a 
normal, unextincted Type la. We describe the main sources of error in such a 
measurement of go, and ways to reduce these errors. 

Subject headings: cosmology: distance scale, dark matter — supernovae: 
general, SN1992bi 
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1. Introduction 

For over 25 years Type I and Type la supernovae (SNe la) have been studied as 
potential "standard candles" for distance measurements (for a review, see Branch & 
Tammann 1992). Recently, Branch & Miller (1993) and Vaughan et al. (1994) emphasized 
the narrow distribution of absolute magnitudes at maximum light for the subset of SNe 
la that have well-measured lightcurves, and are not unusually red or spectroscopically 
peculiar. This subset of 27 "normal" SN la has a dispersion, u Mb = 0.3, that is completely 
accounted for by measurement errors, so the intrinsic dispersion should be smaller still. 
Sandage & Tammann (1993) estimated < j l ^ msic < 0.2 by simulating the effects of intrinsic 
dispersion on Malmquist bias (using a slightly different set of SNe). 

Phillips (1993) noted a correlation between absolute magnitude and lightcurve decay 
time in 9 well-studied SNe la (including a few "peculiars"). Magnitude corrections based on 
such a correlations can "sharpen up" the standard candle by calibrating it, thus improving 
on the current relatively narrow dispersion and also making it possible to include some 
peculiar SNe. For this paper, however, we will take the intrinsic dispersion without this 
calibration to be a l ^ msic ~ 0.25, and we will take the error in the mean absolute magnitude 
for the 27 good SNe la to be <J^ n = a MB /V^7 « 0.06 (note that these are somewhat more 
conservative estimates that those of Sandage & Tammann (1993)). 

Type la's are on average the brightest SNe and therefore could be used to measure large 
cosmological distances and, in particular, the deceleration parameter q of the expanding 
universe (Tammann 1979; Colgate 1979). With many observable features in their spectra 
and their lightcurves, SNe la have an advantage as standard or calibrated candles because 
they can be checked one-by-one for evolutionary differences at high redshift. For example, 
the lightcurve decay time or the blueshift of spectral features may indicate explosion 
strength, and can be checked for systematic differences from nearby supernovae. 
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To use SNe as tools for cosmology, we have developed a strategy to find and study 
high-redshift SNe systematically, using new wide-field, high-resolution CCD cameras and 
image analysis techniques that provide the rapid response time necessary to follow these 
short-lived events (Perlmutter et al. 1993). In 1992 March and April, we demonstrated this 
strategy using the Isaac Newton 2.5 meter Telescope (INT) on La Palma, and found a SN 
at z = 0.458. The farthest previously detected SN was found at z = 0.31 in the SN search 
of Norgaard-Nielsen et al. (1989). 



2. Observations and Data Reduction 

During 1992 March 24 - 28, we observed 54 high-galactic latitude fields, 43 of them 
centered on high-redshift clusters, using 10-minute R band exposures in seeing better than 2 
arcsec. Approximately 200 galaxies were visible in the redshift range z = 0.25 - 0.5 in each 
10.5 arcmin x 11.5 arcmin image to a limiting magnitude of R ~ 23 (4cx detection). During 
1992 April 21 - May 2 the same fields were re-observed, in pairs of 7-minute exposures 
offset by ~6 arcsec to facilitate rejection of cosmic rays and CCD defects. Following each 
night's observations we used semi-automated image analysis software to search for new 
point-sources that had appeared since the reference images were taken. 

With good seeing and 0.57 arcsec pixels, we were able to distinguish candidate SNe 
from other variable sources (e.g. QSO's) by looking for a resolved host galaxy. The time 
between observations at the INT and the completion of analysis in Berkeley was typically 
less than 2 days, so interesting candidates could trigger follow-up observations. The best 
of these SN candidates (SN 1992bi) was found at R.A. = 16 fe 8 m 28!4, Dec. = +39°54'58" 
(equinox 1950.0), 1'.'5 east and 0'.'5 north from the core of an R « 21.3 magnitude galaxy 
(Pennypacker, et al. 1992). We observed the candidate on 5 more nights during the 
2.5 weeks after discovery and on 4 nights during the 9 months after it had faded below 
detectability. Figure [L] (Plate L#) shows the images of the host galaxy before, during, 
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and after the event. To the right of each image is the same image after subtracting off a 
"reference" image from Day 103 of the host galaxy alone. 

The subtraction images of Figure [I] were used to measure the photon flux /|. 2bl of the 
SN candidate on each observation. The Day 396 flux is taken as the best estimate of the host 
galaxy light and is subtracted out of each /|. 2bl measurement. The night-to-night "relative" 
measurement error is approximately 9% at maximum light, and is primarily due to photon 
noise of the background sky. A smaller error contribution, ~5% at maximum light, is due 
to slight mismatches of the two images before subtraction of the host galaxy reference 
image from the measured image of the galaxy plus SN. These errors were determined by 
extensive testing of the plausible range of matching parameters (position offset, seeing, and 
transmission ratio) for their effect on the resulting magnitude measurements. To obtain 
the absolute measurement error we must add (in quadrature with these relative errors) an 
overall error due to the sky noise of the Day 396 reference image, amounting to ~13% near 
peak. Table p] lists for each observation the flux /j| 2bl and relative error of the SN candidate 
measured from the subtraction images (note that /Jf bl can be negative due to photon noise 
of the subtracted host galaxy and subtracted sky). 

The calibration to R magnitude was made on day 149 using the INT/EEV5 camera to 
image M92 and the SN field consecutively (repeated 3 times to check photometric stability) 
with the same airmass within 0.02. These observations were repeated in the V and / bands. 
Our relative photometry agreed with the calibration of 23 stars in M92 (Christian et al. 
1985, revised L. Davis, private communication) with a root-mean-square error of 0.005 mag, 
and yielded the following calibration: m^ 2bl = —2.5 log ffj* hl + 28.20, where /^ 2bl is measured 
in units normalized to the photoelectrons/minute observed at the INT 2. 5- meter on Day 
103. The uncertainty in transferring this standard to each individual night's observation is 
&cai — 0.02 mag, primarily due to the photon noise on ~20 bright stars in the SN field found 
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in both the standard image and each night's images. (Note that the first order extinction 
is calibrated out in this procedure, and that the extinction color dependence is negligible, 
because the calibration observations were in the R band and at an airmass within 0.2 of 
each of the SN observations. The color term is < 0.01 implying that our instrumental r 
band is quite close to the calibration's Cousins R.) 

Twelve attempts to obtain the SN spectrum at four observatories around the world 
were all unsuccessful due to weather (and one instrument failure). On 1992 August 29, 
we obtained two 1800-second spectra of the host galaxy using the Low Dispersion Survey 
Spectrograph 2 on the William Herschel 4.2-meter telescope, and determined the redshift 
z = 0.458 ± 0.001. Both spectra show a strong emission line at 5437 Angstroms. Identifying 
this line as [OII]A3727 at z — 0.458 allows the weak absorption features longwards of it to 
be identified with Ca H, Ca K, and B.S (see Figure 0). 

3. Analysis 

Traditionally one compares photometry of high-redshift objects to photometry of 
nearby objects by calculating a i^-correction to account for the different parts of the 
spectrum that fall in one given filter band for objects at different redshifts. A more robust 
procedure in this case is to compare the photometry of Table [I] to the apparent i?-band 
lightcurve, m^(t), calculated for z = 0.458 from the standard lightcurve of well-observed 
nearby SN la in the B-band. This standard lightcurve, M#(t), is an apparent magnitude 
versus redshift measurement and thus implicitly depends on the Hubble constant chosen. 
However, the expression [Ms(t) — 51ogi^o] that will appear bracketed together throughout 
this Letter is independent of the choice of Hubble constant, as is the resulting measurement 
of go- Since events seen at z = 0.458 last 1.458 times as long as they do locally, this 
lightcurve term must be time dilated to [Mb(£/(1 + z)) — 5 log if q]. The calculated distant 
i?-band lightcurve (for A = 0) is thus: 
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m R (t) = [M B (t/(l + z))-5\ogH }+Am RB + A R + 25 

+ 5 log (cq 2 [l - g + Qoz + (go - l)(2q z + 1)*]) , 

where A R ~ 0.006 is the R extinction in our Galaxy at (l,b) = (63.27,47.24) (Burstein 
& Heiles 1982; Burstein, private communication). The Am R B term is the analog of the 
traditional K-correction: 

Am RB - -^[j FrwSR{x)dX j F SN {X) s B{ X )dX J+2.51og(l + z) 

n. ]n J SFx 0V WS B Wd\ JF™(\')S R (\'(l + z))d\' \ 
g \fFA*V(\)S R (\)d\ JF™(X)S B (X)dX )> 

where S B (X) and S R (X) are the response functions in the B and R bands, and F^ N (X) and 
F^ 0V (X) are the flux per unit of wavelength of the supernova and the standard AOV star 
SED for which B = R = in the Landolt system. 

The correction Am RB accounts for the different zero points of the R and B magnitude 
systems, and for the the redshifting of the SN spectrum. We calculate Am RB = — 0.7 ± 0.05 
using numerical integration with spectra of 3 SNe la with overlapping epochs from 15 
days before to 15 days past maximum (in the SN rest frame). The error represents the 
variation for a given SN over time; the SN-to-SN variation at maximum is 0.02 mag. The 
relatively small variation of Am RB over time is due to the approximate match at z — 0.458 
between the observer's R response function and the SN rest-frame's B response function, 
i.e. Sr(X(1 + z)) S B (X). At z — 0.458, Am RB is therefore less sensitive to the precise 
knowledge of the changing SN spectrum, an advantage over the usual ^-correction. (The 
calculation of Am RB was checked by reproducing the i^-corrections of Hamuy et al. (1993) 
to within 0.001 mag for the 3 SNe.) 

Taking [M B (t max ) — 5 log(i?o/75)] = —18.86 ± 0.06, where the error is the error in 
the mean of the distribution as discussed above (Branch and Miller 1993, Vaughan et al. 
1994), and the template standard B lightcurve, M B (t) — M B (t max ), of Leibundgut (1991 
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and private communication), we calculate the lightcurves shown in Figure FJ for two values 
of the deceleration parameter, and 0.5. At z = 0.458 a standard candle is 0.25 mag fainter 
(for A = 0) in an empty universe, go = 0, than in a critically closed universe, go = 0.5. 

Figure ^ also shows the measured photometry data from Table [j]. When we perform 
a x 2 fit °f th e th e data points of Table [I] to the calculated lightcurve, fn(t), varying the 
values of go and the time offset, At, we find the best fit at go = 0.07. If there is significant 
host galaxy extinction then g would be larger. The individual error bars shown in Figure fj 
are the relative errors used in the fit, while the errors from the reference image photometry 
and calibration contribute to the overall error bar o c (discussed below). Note that the pre- 
and post-SN measurements constrain the time-offset fit. 

Although the photometry data shown in Figure fj are clearly consistent with the 
lightcurve of an SN la at z — 0.458, several alternative identifications were considered. It 
is very unlikely (Prob < 1%) that this event is a projected foreground variable star in our 
Galaxy or a projected background quasar or AGN. The lack of repeated flare events makes 
these identifications still more unlikely. The data are also consistent with an unusually 
bright SN IIL. Assuming the currently favored cosmologies with go < 2 (or, equivalently, 
H > 40 km sec -1 Mpc -1 and r un i verse > 12 Gyr; see Carroll, Press, & Turner 1992), we 
can estimate the relative probability of a SN IIL identification as follows: First, we note 
that the event's absolute magnitude could not be more than 0.7 magnitudes fainter than 
the SN la value M B {t max ) used above, because g would then have to be > 2 to fit our 
photometry data. We then compare, in nearby spiral galaxies, the number of SNe la and 
the number of SNe IIL brighter than this minimum: M m i n = Msi^max) + 0.7 + 0.3 (the 
additional 0.3 magnitude allows for measurement error). The result is ~10 SNe la for 
every SN IIL (see Miller & Branch 1990), and we use this ratio as the relative probability 
of these identifications. This is, of course, a rough approximation that assumes that the 
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SN detection efficiency falls off towards fainter absolute magnitudes at least as fast at 
z = 0.458 as in nearby searches. Similarly, we estimate < 1/30 for the relative probability 
that this event is a peculiar SN la or one of the SNe la showing clear evidence of extinction, 
since over 30 times as many normal, unextinguished SNe la brighter than M min are found 
in nearby galaxies (Branch & Miller 1993). The identification of this event as a normal 
SN la at z — 0.458 is therefore made with ~90% confidence level (for any q < 2) for the 
purposes of this paper. The confidence level would be greater (and the assumption of go < 2 
unnecessary) if spectra or color photometry were available. 

The uncertainty in this measurement of go has two major contributory sources. The 
first source is the uncertainty, cr R ~ 0.14, in the apparent R magnitude at maximum light 
of the z = 0.458 SN. This includes the photon noise and image-matching uncertainty in the 
images observed near maximum light, cr peak ps 0.06, and the photon noise and calibration 
uncertainty in the Day 396 reference image used to subtract off the host galaxy light, 
cr re f ps 0.13. (A better reference image will reduce this error source significantly.) 

The second and much larger contributory source of measurement error enters in when 
this distant photometry is compared with the nearby apparent magnitude-Hubble velocity 
relation, i.e. the nearby measurement of [Mb — 51ogi?o]- The uncertainty in the mean of 
this quantity may be as small as cr](J? an ~ 0.06, but since we do not know where in the 
M B distribution our particular SN at z — 0.458 lies, we must add in quadrature an error 
contribution for the intrinsic dispersion a l ^ insic ps 0.25. This is currently the dominant 
source of error. Since it is plausible that this is an overestimate of the intrinsic dispersion, 
as discussed above, we will separate out this error in the subsequent analysis, and define a 
a Q to include only the other sources of error. Finally, we also include in quadrature an error 
of a RB ps 0.05 for the uncertainty in Am RB , so — a R + <7]JJ? an2 + a RB . 

The total error in the distant SN apparent magnitude and in the comparison with the 
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nearby apparent magnitude-velocity relation is a mag = ±a ± < j l ^ ms%c = ±0.16 ± 0.25. 
(Figure |3|'s a c is the equivalent total error in flux.) This propagates through to an uncertainty 
on the go measurement of a qo ~ <T mag /z, yielding a measurement of go = 0.1 ± 0.3 ± 0.55 
(for A = 0). This measurement is stricter in its lower limit than its upper limit, since host 
galaxy extinction would increase go- We emphasize again the assumption of a normal SN la. 

4. Discussion 

We wish to stress two main results: 

(1) From a single SN at z = 0.458 known to be a "normal" type la, we could in 
principle measure go with an uncertainty comparable to the previous measurements in 
the literature that required larger numbers of objects (galaxies or clusters) and needed 
significant corrections for evolutionary effects (for reviews, see Rowan-Robinson 1985 and 
Sandage 1988). Clearly the measurement of go cannot rest on a single distant SN la, and in 
this particular case we do not have the color photometry or the spectra that would enable 
us to screen for host galaxy extinction and peculiar SNe. We are continuing to search for 
distant SNe la, and we are scaling up the detector size and sensitivity so that we can find 
many high-redshift SNe la, and follow their lightcurves, colors and spectra over maximum 
light. This will make it possible to compare the luminosity function distribution of distant 
and nearby SNe, and to identify and reject extinguished SNe. 

(2) The measurement errors that contribute to the uncertainty in g can all be reduced. 
The dominant error is due to both our uncertainty in the true width of the SN la absolute 
magnitude distribution and our uncertainty in where a given distant SN la falls in this 
distribution. As a few more SN la past Coma are discovered and studied, we may find 
that the distribution width is smaller than can be seen in the nearby SNe where peculiar 
velocities can mask the Hubble expansion. Hamuy et al. (1994) are working on such a 
search, and the Berkeley Automated Supernova Search (Muller et ai, 1992), now being 
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moved to a good astronomical site, may also find and study these SNe. 

Alternatively, if the dispersion is in fact <j^* rmsjc = 0.25 then discovering N more 
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distant SNe la will make it possible to statistically reduce this error by ^V^/V- The main 
source of concern here is that the distant M B distribution should be the same as the 
nearby distribution for the comparison to be valid; it will be important to ensure that the 
distant SN detection limit be significantly fainter than the peak of the M B distribution at 
the redshifts searched to avoid Malmquist-bias distortion of the shape of the distribution, 
particularly if there are rare superluminous Type la's (see the discussion of the possibly 
superluminous SN 1991T in Ford et al. 1993). 

The uncertainty in the .fT-correct ion (in this case, Am^) is the only remaining source 
of error that cannot be improved simply by using longer exposure times, larger telescopes, 
or darker sites. It therefore may be useful to measure the lightcurves for future nearby SNe 
la using filter sets designed to give "blue-shifted R" bands for a sampling of values of z. 
Together with SNe la spectra to interpolate between these z values, these "blue-shifted R" 
lightcurves should allow accurate determinations of the X-corrections. 

It will be necessary to check the effects of metallicity and other evolutionary changes 
on the SN la absolute magnitude to validate measurements of q derived from SNe la. 
Deep spectra of distant SNe will be necessary to look for any evolutionary changes with 
respect to spectra of nearby SNe. Metallicity effects can be studied by comparing nearby 
SNe la found in different metallicity environments. So far, however, no indications of 
magnitude evolution have been found observationally, and the standard model of SNe la 
gives theoretical reasons to believe that SN la magnitudes should not evolve. Thus, the 
prospects for obtaining a precise value for q from SNe la appear good. 

[Note added in press: This project has by now found 6 additional SNe around z ~ 0.4 
(e.g., Perlmutter et al 1994). With these it should be possible to use lightcurve decay-time 
calibration or lightcurve shape calibration, recently shown to yield a M t " nsic as small as 0.1 
(Hamuy et al. 1994) to 0.2 (Riess et al. 1994).] 
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Day 
JD-2448707 


Telescope 
/ Camera 


Exposure 
(seconds) 


x92bi 

Jr 

(relative fTux) a 





INT/EEV5 


600 


-18 ±42 


2 


INT/EE V5 


600 


25 ±25 


30 


INT/EE V5 


2x420 


260 ± 21 


34 


INT/EE V5 


2x420 


251 ±22 


38 


INT/EEV5 


2x480 


241 ± 29 


43 


INT/EE V8 


2x300 


213 ±23 


45 


INT/EE V8 


300 


169 ±41 


47 


INT/EE V8 


300 


215 ± 34 


77 


Palomar/COSMIC 


300 


28 ±44 


103 


INT/EEV5 


5x~170 


20 ± 19 


148 


INT/EEV5 


6x300 


12 ±22 


330 


INT/Ford 


6x320 


-34 ± 22 


396 


INT/EEV5 


7x400 





Table 1: Observation log and supernova flux 



a The host galaxy flux from the image of Day 396 has been subtracted off. Relative flux units 
are normalized to the photoelectrons/minute observed at the INT 2.5-meter on Day 103. 
The calibration is: m 9 ^ hl = — 2.5 log /]f bl + (28.20 ± 0.02). The listed relative error includes 
sky photon noise and image matching error for each individual image. An overall error of 30 
units should be added in quadrature to account for the sky noise of the reference image. 
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FIGURE CAPTIONS 



Figure 1 (Plate L#): The left panel of each pair shows the image of the host galaxy 
with or without the supernova. The time, t, is JD — 2448707. The superposed contours 
start at 22 photoelectrons/minute/pixel, with each additional contour representing 11 
photoelectrons/minute/pixel. Note that these images are normalized to the same sky 
transmission but not matched in seeing. The right panel of each pair shows the same 
images after subtracting off the host galaxy image of day 103. (Days 0, 330, and 396 are 
not shown.) 

Figure 2: Spectrum of the host galaxy, observed on 1992 August 29 (Day 156) on the 
William Herschel 4-meter telescope. 

Figure 3: Solid curve shows the calculated lightcurve, fii(t), for the best fit g = 0.1, 
based on the template B lightcurve for nearby supernovae. The dotted curves are fn(t) 
for go = 0.5 (upper curve) and for go = (lower curve). The photometry points are /|j 2bl 
from Table |l| (days 330 and 396 not shown). The inner error bar, cr c , shows the combined 
uncertainty at maximum light in fn(t) and in the reference image's photometry and 
magnitude calibration; the outer error bar includes the intrinsic dispersion, u l ^ msic _ The 
inset scale shows R magnitudes. For the best fit go, the peak magnitude is ttir = 22.2 on 
Day 34. 
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